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Abstract 

We study the quantum Bethe ansatz equations in the 0{2n) sigma- model for physical par- 
ticles on a circle, with the interaction given by the Zamolodchikovs' S-matrix, in view of its 
application to quantization of the string on the 5^""^ x Rt space. For a finite number of 
particles, the system looks like an inhomogeneous integrable 0(2ra) spin chain. Similarly to 
OSp{2m + n\2m) conformal sigma-model considered by Mann and Polchinski, we reproduce 
in the limit of large density of particles the finite gap Kazakov-Marshakov-Minahan-Zarembo 
solution for the classical string and its generalization to the x Rt sector of the Green-Schwarz- 
Metsaev-Tseytlin superstring. We also reproduce some quantum effects: the BMN limit and 
the quantum homogeneous spin chain similar to the one describing the bosonic sector of the 
one-loop = 4 super Yang-Mills theory. We discuss the prospects of generalization of these 
Bethe equations to the full superstring sigma-model. 
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1 Introduction 



Superstring theory on non-trivial gravitational backgrounds is extremely important for the 
progress of the theory of fundamental interactions and cosmology: not only it might describe 
some rather realistic models in cosmology and black hole physics but it also provides a dual 
description of some strongly interacting planar Yang-Mills gauge theories. A set of these du- 
alities are often named as AdS/CFT correspondence [1-3] (see [4] for the review). Apart from 
its physical significance for the strongly interacting gauge theories a great technical advantage 
of the string side of duality is that the string theory in the tree approximation is a two dimen- 
sional (T-model, at least in the world-sheet perturbation theory, and the string interactions are 
not relevant in the planar 't Hooft limit of the dual gauge theory. And since the two-dimensional 
cj-models are sometimes solvable it gives us hopes that so are some realistic string sigma-models, 
together with the dual supersymmetric Yang-Mills (SYM) theories in the large N limit. 

The master-example of the AdS/CFT correspondence, the 4-dimensional conformal J\f = 4 
SYM theory, and its dual, the Green-Schwarz-Metsaev-Tseytlin superstring [5] on the AdS^ x 
background looks particularly promising in this respect. The 1-loop integrability was discovered 
in = 4 SYM in [6] for the bosonic sector^ where the dilatation operator was identified 
with the Hamiltonian of an integrable 1-dimensional spin chain. It was soon followed by the 
demonstration of the classical integrability of the full superstring cr- model on AdS^ x [9], 
where already some exact particular classical string solutions were known [10-12]. It was then 
shown in [13] that the integrability might persist in = 4 SYM even in higher loops, and is 
certainly a property of the full theory at one loop [14, 15]. The underlying (super) spin chains 
are solvable by the standard Bethe ansatz technique and possess an interesting ferromagnetic 
continuous limit [16-18] corresponding to very long local single trace operators in SYM. The 
resulting anomalous dimensions were directly compared for some states with the energies of 
classical string solutions in the so called Frolov-Tseytlin limit. 

As a next step in this progress, the general classical solution for the superstring rotating on 
AdS5 X was proposed, first for the 5^ x Rt sector [19], then for AdSs x [20], x Rt [21], 
and finally for the full AdS^ x superstring [22]. The full solution even bears an imprint 
of fermions in the corresponding algebraic curve describing a finite gap solution. It was shown 
in [19,21,23,24] that a// states and anomalous dimensions of the SYM theory can be mapped onto 
all classical states of the string in this particular limit, giving another impressive demonstration 
of the AdS/CFT correspondence. 

This integrability is recently generalized to three-loops [14, 25] and even a plausible pertur- 

bative all-loop ansatz was proposed [26,27]. The perturbative comparison of the AdS/CFT 

correspondence works well for this ansatz up to two loops but breaks down at three loops [28], 

a mismatch usually explained by certain non-perturbative effects which are difficult to take into 

account in the comparison of the weak and strong coupling regimes.^ One- loop world-sheet 

^Integrable spin chains have been discovered in (non-supersymmetric) gauge theories before [7,8]. 
■^Nevertheless, the classical string KMMZ solution [19] inspired a beautiful result of [25], where the three-loop 
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calculations for particular string solutions [29], as well as the direct one-loop world-sheet calcu- 
lation of the BMN spectra [30-32] also show a perfect AdS/CFT correspondence with respect 
to similar quantum corrections to the solutions in SYM spin chains in the scaling limit [33-36]. 

In spite of the multitude of signs of integrability of the = 4 SYM, the ultimate progress 
is hardly possible without the establishment of integrability of the full quantum string cr-model. 
It is not an easy task, though far from being hopeless (see a discussion in [37]). Similar a- 
models were already proven to be integrable starting from the Zamolodchikovs' discovery of the 
factorized unitary and crossing-invariant S-matrix for physical particles in the 0{n) sigma- model 
[38]. Then the SU{n)L x SU{n)ji principal chiral field was solved [39-44] and the corresponding 
physical S-matrix was found [40] . This model represents the first non-trivial example of the field 
theory completely solvable in the 't Hooft limit [45]. Some supersymmetric o"- models, including 
conformal ones, appeared to be integrable as well [46-48]. 

The main difficulties for establishing the integrability of the full superstring are the following: 

1. It is a non-compact u- model as the AdS^ subspace is. The S-matrix approach is compli- 
cated by the existence of a continuous spectrum of asymptotic states, although we know some 
integrable cr-models of this kind: quantum Liouville theory [49-52] bearing many similarities 
with a coset of AdS^ sigma-model [53], the Sine-Liouville theory and its strong-coupling dual, 
the 2D dilatonic black hole ^. 

2. The Virasoro conditions impose a complicated selection rule on the possible quantum 
states: they eliminate the longitudinal motions of the string forbidden by the world-sheet 
reparametrization symmetry, together with the corresponding quantum states. We could start 
from a gauge which eliminates the conformal redundancy from the beginning, but it hides the 
world-sheet Lorentz covariance as well. 

3. It seems hopeless to find a completely closed subsector of the superstring on AdS^ x at 
the quantum mechanical and nonperturbative level, and to write down a selfconsistent system 
of Bethe equations for it (though see some instructive attempts in [58-60] ) . To obtain the right 
divergency free string sigma model we have to consider the full theory, including fermions. 

The full Metsaev-Tseytlin superstring looks like two sigma models, one on AdS^, another 

on space, each containing 4 on-shell bosonic fields "glued" together by 16 fermionic fields (8 

on-shell fields) . Classically the fermions do not influence the dynamics and the system looks like 

two independent bosonic cr-models. If we quantize the cr-model (with a trivial inclusion of the 

AdS-time coordinate, which makes it x Rf) it will show the asymptotically free behavior, even 

if we impose the analogue of Virasoro constraints (no longitudinal modes). Nevertheless, this 

compact 0(6) cr-model may provide a valuable information on the full superstring: its classical 

limit describes precisely the compact sector of the superstring [21] and its quantization should 

bear some important features of the full quantum superstring cr-model. And if we manage to 

dilatation operator was constructed for some subsectors of N=4 SYM. 

■^Similarly to the AdS/CFT correspondence of superstring to the SYM, they also have their dual, matrix 
models and matrix quantum mechanics description for the c = 1 matrix model [54], or the matrix black hole 
model of [55] (see [56] for the discussion of this issue). 
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impose the Virasoro constraints in a natural way it can serve as a rather sophisticated toy model 
for the quantum superstring. 

Our principal motivation in this paper is to study the full quantum version of the 0(2n) 
sigma model in the context of its application to the string theory. Our starting point is the 
Zamolodchikovs' S-matrix for physical particles. We put L such particles on a circle and study 
possible states by means of the nested Bethe ansatz equations diagonalizing the underlying 
monodromy matrix by the use of the quantum inverse scattering method [61,62]. The main 
result may be interesting not only to the string theorists but also to the specialists in integrable 
2D field theories: we reconstruct all finite gap solutions of [19,21] for the classical string motions 
on X Rt, directly from these Bethe equations in the limit of large density of the particles with 
sufficiently big momenta. In this limit, the theory is in the extreme asymptotically free regime. 
We can neglect the conformal anomaly since the particles are almost massless, and the states 
are classical solitons. These solitons are nothing but the classical string states. The string looks 
like a collective state of the particles. The string mode numbers can be identified with certain 
"magnon" mode numbers of the Bethe equations, the number of excited modes corresponds to 
the number of gaps in finite gap solution (this number is equal to the genus of the algebraic 
curve describing the solution), and the amplitude of a mode is fixed by the number of Bethe 
roots forming the corresponding gap. 

The system looks like an inhomogeneous 5*0 (6) spin chain with the nearest neighbor inter- 
actions defined by the rapidities of the neighboring particles. The rapidities, in their turn, are 
fixed from the periodicity condition, thus make the spin chain dynamical^. 

Qualitatively it is quite similar to the OSp{2m + n\2m) supersymmetric u-model of [64], 
conformal in the limit n — > 2. Our way to restore the classical limit is also inspired by [64], 
but the details are very different, since our model is conformal only classically. However, some 
features of our model look closer to compact sectors of the string theory on AdS^ x S^. In 
particular, in a certain limit we restore the homogeneous spin chain, the same (up to some 
normalizations) as the one describing the 1-loop SYM theory in the corresponding sector. To 
make it working for higher SYM loops it might be necessary to incorporate into the future string 
(T-model certain features of the Hubbard model which is a good candidate for reproducing all 
loops of the perturbative SYM [65]. 

The classical string solutions are encoded into an algebraic curve. To each curve with the 
fixed moduli one can associate the unique (up to some trivial choice of initial conditions, as 
proven in [66]) classical string solution with Virasoro conditions imposed and with the number 
of excited oscillatory modes corresponding to the genus of the curve. 

The algebraic curve obtained in [19,21] from the purely classical finite gap method [67] turns 

out to be the same as the one extracted here from the Bethe equations, and the two different 

''Not to confuse it with the dynamical S-matrix of [63] for the su(3|2) sector of the J\f — 4 SYM where even 
the length of the chain can change, a feature certainly to be taken into account for the full string cr-model. 
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projections (Riemann surfaces) of the curve are related by Zhukovsky map 

Z = X + 1/x ^ 

which is an important part of the construction^. It was used aheady in [19] for the successful 
2-loop comparison of the string and SYM results and its significance is better understood in 
[26,58,66]. This map might identify the right variables to quantize in the full superstring theory. 

The Virasoro conditions deserve a special mentioning here. As it will be seen from our 
classical limit procedure, the correct classical string solutions are related only to Bethe states 
where all rapidities of the physical particles have the same mode number. In the classical limit 
they condense into a single cut, which reproduces, after Zhukovsky transformation, the two pole 
structure of the algebraic curve of the finite gap solution. We conjecture that this is the way 
to impose the Virasoro-like constraints for any general quantum state as well: the rapidities 
should have the same mode number. We confirm it by some examples of particular solutions 
with small longitudinal amplitudes.^ 

Finally, we would like to mention that the Bethe ansatz description provides an unusual 
example of the holography in the string theory: the Bethe roots look like the positions of some 
D-branes reducing the quantization procedure to a certain "discretization" of the underlying 
algebraic curves. But unlike the matrix model description of D-branes where the eigenvalues 
serve as quantum coordinates of the D-branes and should be dynamical variables, the Bethe 
roots of the stringy o"-models are fixed by the periodicity condition. 

The paper is organized as follows. In section 2, we define the 0{2n) cj-model and remind 
its classical integrability properties allowing the finite gap solution. In section 3, the quan- 
tum Bethe ansatz solution of the SU{2) principal chiral field model, equivalent to the 0(4) 
cj-model, is presented and its classical limit is considered and the corresponding algebraic curve 
is constructed. In section 4, we map this algebraic curve to the curve of the finite gap solution 
obtained directly from the classical action, fixing appropriately the relations between all the 
parameters of two approaches. In section 5, the construction is generalized to the 0(6) sector. 
In conclusions we make some remarks concerning possible generalizations of this construction 
to the full superstring on AdS^ x and problems related to it. In Appendix A the solution 
for the U{1) sector in the scaling limit is presented in detail. In Appendix B we provide an 
alternative proof of our correspondence for the 0(4) sector, namely through the transformation 
of resolvents. In Appendix C both the bootstrap method as well as the algebraic Bethe ansatz 

^Thc Zhukovsky map was proposed in 1906 in [68] for description of aerodynamics of the airplane wing. It 
maps a circle passing through the point x — 1 into a wing-shaped figure in z-plane (see the front page). 

^At the classical level Virasoro constraints have a clear interpretation as a restriction on the allowed solutions 
forbidding longitudinal motions of the string, even if we truncate the full AdS^ x theory to a subsector, like 
the X R subsector considered in this paper. On the quantum level, the subsectors are not even described by 
a conformal field theory and the Virasoro constraints do not have a clear meaning and are not imposed by the 
reparametrization symmetry. Therefore we rather have to speak in this case about a natural selection rule for 
the quantum states which leads to the standard Virasoro conditions in the classical limit. 
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program are reviewed for the 0(4) o"- model. In Appendix D this is generahzed to the 0{2n) 
fj-model. 



2 0(2n) Sigma-Model and Classically Integrable String Theory 

The main system of our interest is the superstring on the AdS^ x background. A compact 
bosonic subsector of it is described by the sigma model on the subspace x Rt, where Rt is 
the coordinate corresponding to the AdS time. The time direction will be almost completely 
decoupled from the dynamics of the rest of the string coordinates, appearing only through the 
Virasoro conditions. These conditions represent only a selection rule upon the states of the 0(6) 
sigma model allowing to choose a particular set of states, or a particular set of solutions in the 
classical limit, which do not contain any dynamics along the string. 

Of course in the absence of fermions and the AdS part of the full lOD superstring theory, 
this model will be asymptotically free and will not be suitable as a viable quantum string theory. 
In addition, on the dual = 4 Yang-Mills side, the corresponding 50 (6) sector of bosons is not 
closed under the action of dilatation operator (except the limit of long operators which can be 
compared with the classical string [69]). However, in the classical limit we will encounter the 
full classical finite gap solution of the string in 5*0(6) sector found in [21]. Furthermore we still 
expect to capture some features of the full quantum string cr-model. 

In the rest of this section we will first take instead of 0(6) a more general 0(2n) case. We will 
later reproduce the full classical finite gap solution [21] of this model by means of the Zamolod- 
chikovs' physical S-matrix which may be considered as the most convincing demonstration of 
the correctness of the quantization of the model by the bootstrap method. 

2.1 Definition of the Sigma-Model on S'^^~^ and Virasoro Conditions 

The X Rf reduction of Green-Schwarz-Metsaev-Tseytlin superstring on AdS^ x background 
[5, 70] in the orthogonal gauge has a simple action [12] in terms of homogeneous target space 
coordinates Xj(r, cr), i = 1, . . . , 2n and a scalar Y{t, a) representing the AdS time: 



The coupling constant in front of the action is identified by the AdS/CFT correspondence with 
the 't Hooft coupling A = g'^Nc of the = 4 supersymmetric Yang-Mills (SYM) theory. 
We work in the static gauge Y{T,a) = kt. The Virasoro condition reads 




(2.1) 



d± = dr± da 



(2.2) 



where the current is a 6 x 6 matrix defined as 



{j±)ab = {h~'d±h)^^ = 2 {Xa d±Xb - d±Xa Xfe) , hab = d^b " 2V,Vfe. 



(2.3) 
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K is identified with the dimension of the corresponding SYM operator according to the 
AdS/CFT correspondence. 

To ihuminate the difference between the string sigma-model on x Rt and the ordinary 
sigma- model on 5^, it is instructive to consider shghtly more general constraints. These con- 
straints can be viewed as the Virasoro constraints for another gauge 

Y='^K+{T + a) + ^K.{T-a). (2.4) 

AdS time Y is periodically compactified in this case. 

tr(j±)' = -2 4. (2.5) 

The complete Virasoro conditions forbidding the time compactification will be imposed by the 
identification = k„ = k. Note that the condition eq.(2.5) still allows a nonzero longitudinal 
momentum, although its density, as well as the energy density are constant along the string. 
In what follows we will deal with the action 



tr {jaja) , (2.6) 



subject to the condition eq.(2.5), with an arbitrary number of fields Xi, i = 1, . . . , 2n. 
The energy and the momentum are given by 

E'''±P''' = -^ld'atr\jo±n?, (2.7) 

or, using (2.5), 

E^^ = ^{.l + .l), P^i = ^(4-.2_). (2.8) 
2.2 Classical Integrability and Finite Gap Solution 

The equations of motion and the form of the current, j = h^^dh, can be encoded into a single 
fiatness condition [67] 

[£+(x),£_(x)] =0, (2.9) 
where the Lax pair of operators, the currents deformed by spectral parameter x, are given by 

C±ix) = d±-^. (2.10) 
The Lax operators describe a connection over the world-sheet and define the monodromy matrix 

exp£'d4(-j=?i^ + ^). (2.11) 

By construction i}{x) is a complex orthogonal matrix and thus has the eigenvalues 

^g«9i(a;)^ g-*gi{a;)^ gi(?2(2:)^ g-ig2(a;)^ _ _ _ ^ giqn{x) ^ ^-iqn{x)^ (2.12) 
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where Qkix) are called quasi-momenta. They do not depend on time r due to eq.(2.9) and 
provide an infinite set of classical integrals of motion of the model. As for the global conserved 
charges of the sigma-model one has 

J=- — : / dcrjV = diag(Ji, -Ji, J2, -J2, . . . , Jn, -Jn) • (2.13) 

iTTl J 


This current is normalized so that after quantization Jj are integers for any quantum state. This 
is the 0(2n) angular momenta quantization condition. 

Let us now recall some analytic properties of these quasi-momenta [21]. Expanding the 
monodromy matrix in powers of 1/x around x = 00, we obtain 

q,ix) = -^ + 0{l/x'). (2.14) 

We see from (2.11) that simple poles will appear at x = ±1 for the quasi-momenta. Furthermore, 
from (2.3) we see that j has only 2 eigenvectors corresponding to nonzero eigenvalues. Thus we 
can diagonalize j± in such a way that the poles at ±1 are only present in qi, 

qkix)=Sk,i^^ + 0{ixTlf) forx^il. (2.15) 

Furthermore the monodromy matrix (2.11) with inverse spectral parameter, ^}{l/x), can be 
written, using the explicit form (2.3) of j and h, as 

n(l/x)=Pexp [ da (h-^- ( + ^—^ h- h-^d^h] = h(2Tr)n(x)h-^ (0) . 
Jo V 2 Va; - 1 x+ ly J 

which is nothing but a similarity transformation due to the periodicity /i(27r) = /i(0). Therefore 
S7(x) and f2(l/x) have the same set of eigenvalues. For the quasi-momenta this gives the 
transformation law 

qk{l/x) = dk,i (47rm - qi{x)) + (1 - Qk{x) (2.16) 

which respects (2.15). This symmetry will be of upmost importance in our further discussions. 
From this relation one immediately reads off the analytic behavior of ijk around x = giving 
the behavior around x = 00. 

Finally we obtain the quasi-momentum by solving the characteristic polynomial of which 
is well defined in C. This, on the other hand, implies that the quasi-momenta will be the 
2n branches of a single function with branch cuts Ca- Along these cuts the quasi-momenta in 
general get permutated and shifted by multiples of 27r as 

h^ii = '^T^na , x£Ca, (2.17) 

so that we obtain some sheets connected by the cuts as in the example in figure 6. For each cut 
outside the unit circle there will be a mirror cut inside the unit circle due to (2.16). 
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2.3 Classical SU{2) Principal Chiral Field 

In this subsection we will concentrate on the classical finite gap solution of the 0(4) sigma 
model formulated in terms of the SU{2) principal chiral field. We will essentially briefly repeat 
the construction of [19] (with a small generalization to the excitations of both left and right 
sectors) to fix the notations for the easy comparison with the quantum Bethe ansatz solution 
of the model. 



2.4 Definition of the Model 

The 0(4) sigma model represents the reduction of the sigma model in the AdS^ x background 
to the subsector of string moving on x Rt. Classically this is a perfectly consistent reduction 
while at the quantum level one might still expect to capture some features of the full theory. 
The action of the theory can be represented in terms of the SU{2) group valued field g = 
Xi +iTsX2 + iT2X3 + iTiXi {n are the Pauli matrices). Then the action eq.(2.6) takes the form 
of the SU (2) principal chiral field 



Svr 



da j tv{dag^dag) . (2.18) 



The obvious global symmetry of the action is the left and right multiplication by SU{2) group 
element. The currents of this symmetry are 

j±=9~'d±g, ji = d±gg-\ (2.19) 

with corresponding Noether charges 

Ql = ^ j da tr (^005 9^=^) , Qr = ^ j da ti (ig^dogr^) . (2.20) 

In the quantum theory these charges are positive integers. It will be important for future 
comparisons to notice that the normalization of the generators is such that the smallest possible 
charge is 1. 

Virasoro conditions read 

tr(i2) = _24. (2.21) 
From the action we read off the energy and momentum as 

E^'±P'^' = -^ /trbo±Jifrf^ = ^4- (2.22) 

The Lax construction (2.9-2.11) can be repeated for this formulation using either the left or 
the right current, the two choices being related by a simple relation 

nUl/x) = g{2TT)nRix)g~H0) . (2.23) 
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Let us use Hr, for definiteness. Now there will be only two quasi-momenta pi{x),p2{x) since 
the 2x2 unitary monodromy matrix will have eigenvalues e*^^^^\ e*^^^^^ . Unimodularity 
imposes 

Pi{x) = -p2{x) mod 2tt (2.24) 

so that we can define 

T{x) = ]^T:vVLr{x) = cosp(x) . (2.25) 
Analyzing the singularities of eq.(2.11) we find the behavior of the quasi-momentum at x — > cxo 

p(.) = -^ + ..., (2.26) 

VAX 



and at X — > (from eq.(2.23)) 



p{x) ~ 27rm + ^^!^x + • • • , (2.27) 
v A 



where we also used the fact that due to the periodicity of (?(ct), 17(0) = 1 and hence p(0) = 27rm. 
Finally, at x — > ±1 

m = -^ + --- , (2.28) 

X =F 1 

By construction 0/j(x) is analytical in the whole plane except x = ±1 where one has essential 
singularities. Then, from eq.(2.25), one concludes that for x ^ ±1 the only singularities of 

can be of the type p' {x — > x^) ~ / ^_ . ■ If we look for the finite gap solution the number of 
these branch points is finite and even, and we conclude from (2.26-2.29) that the Pi{x),p'2{x) 
are two branches of an analytical function defined by a hyperelliptic curve of genus g: 

The coefficients of the polynomial Pg-i (x) of degree g — 1 and b±,c±, and the positions of branch 
points Xfc make a total of 3g' + 6 constants. They are fixed by 

1. The asymptotics (2.26-2.28) and the absence of simple poles at x = ±1 (5 conditions). 

2. The singlevaluedness conditions 

dp = j = l,...,ff (2.31) 



where the homomorphic integrals around the cycles Ai are essentially the contour integrals 
around the cuts Ck {g conditions). 
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3. The integer S-period conditions, analogous to eq.(2.17) 

dp = 27mj j = l,...,g + l. (2.32) 



following from the ambiguity of the definition of the branches of the quasi-momentum 
eq.(2.24) {g + I conditions). 

The cycle Ai is the contour which encircles the cut Ci - see figure 4 - in the counterclockwise 
direction in the upper sheet if the cut is outside the unit circle or in the clockwise direction in 
the lower sheet if it is inside the unit circle. 

The -Bfc cycle consists of the contour starting at 00+ in the upper sheet, changing sheet in 
at the cut and going from here to oo_ in this lower sheet. 

We can also write (2.27) as 

dp = 27rm . (2.33) 



I Be 

where this new cycle goes from 00 to in the upper sheet. 

Let us distinguish the cuts inside and outside the unit circle by superscripts u and v. The 
remaining g constants can be parametrized by the filling fraction numbers which we define as 

= £. P(-) (1 - ^) Sr = ^ (1 - J,) (2.34) 

From the AdS/CFT correspondence these filling fractions are expected to be integers since this 
is obvious on the SYM side [19,24]. In fact, it was pointed out in [21,24] and shown in [66] 
that S*"'^ are action variables so that quasi-classically they indeed become integers. We also 
find a striking evidence for this quantization in the string side when finding the classics from 
the quantum Bethe ansatz where these quantities are naturally quantized. 
The conditions eq.(2.33) can also be represented as 

#(x) = 2^nfc, x€Ck (2.35) 

where we defined 

Hx) = \{p{x + iO)+p{x-iO)). (2.36) 

3 Quantum Bethe Ansatz and Classical Limit: 0(4) Sigma-Model 

We will describe a quantum state of the 0{n) sigma model (the analogue of the "closed string" 
state for the x Ri sector of the superstring theory) by a system of L physical particles of the 
mass TTiQ put on a circle of the length C These particles transforms in the vector representation 
under 0(4) symmetry group or in the bi-fundamental representations of SU{2)ji x SU{2)l. The 
scattering of the particles in this theory is known to be elastic and factorizable: the relativistic 
S-matrix S {61 — O2) depends only on the difference of rapidities of scattering particles 61 and 62 
and obeys the Yang-Baxter equations. As was shown in [38] (and in [40,41,43,44] for the general 
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principle chiral field) these properties, together with the unitarity and crossing-invariance, define 
essentially unambiguously the S-matrix. In Appendix C we review the Bootstrap program of [38] 
as well and the algebraic Bethe ansatz construction for this system [71]. 

There should be no confusion with the employed notation. Having in mind the comparison 
in section 4 we will use the same letters in sections 2 and 3 for quantities that we be later proven 
to be the same. For what concerns this section they are just definitions. 

3.1 Bethe Equations for Particles on a Circle 

When this system of particles is put into a finite 1-dimensional periodic box of the length L the 
set of rapidities of the particles {^a} is constrained by the condition of periodicity of the wave 
function \'\\)) of the system, 

1 N 

where the first term is due to the free phase of the particle and the second is the product of 
the scattering phases with the other particles. Arrows stand for ordering of the terms in the 
product. [1 = moC is a dimensionless parameter. The physical SU (2) PCF S-matrix is a product 
of left and right S-matrices, S = Sl x Sr [38, 40] 

and 

It has the following large 9 asymptotics: 

iiogsi{e)^i/e + o{i/e^). (3.3) 

The diagonalization of both the L and R factors in the process of fixing the correct periodicity 
(3.1) leads to the following set of Bethe equations 



g— i/x sinh irBa 



TT Q2/fl a ^ TT - + IT " + (o a\ 

_ -rr Uj -9p- i/2 ^ Uj-Uj + i 

- llu,-9, + ^/2l\u,-u.-^^ ^'''^ 
1 = U '^'-l'-J' U'-^^^^ (3-6) 

where 

a, /3 = 1, . . . , L, i, j = 1, . . . , Ju, k,l = 1, . . . , Jy . 
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The left and right charges of the Bethe vector, i.e. wave function, associated with the two SU{2) 
spins are given by 

Ql = L — 2Ju , Qr = L — 2 Jt, , (3.7) 

since a u or v root correspond to a spin flip in the corresponding SU (2) and since we are 
normalizing charges so that 1 is the unit charge. The absence of n's and v^s corresponds to the 
ferromagnetic limit with all spins aligned in the up direction. 

Both the construction of the ^-matrix by means of the bootstrap program [38] and these 
Bethe equations (see a similar example in [71]) describing the solution to (3.1) are reviewed in 
Appendix C for the sake of completeness. 

The total momentum and the total energy of the state are the sums of momenta and energies 
of individual particles. In accordance with eq.(2.18) we take £ = 27r (in the units of the string 
tension). Then 

^' = ^^sinh(7r0,), (3.8) 

^ = 77^Vcosh(^^,) (3.9) 

This model with massive relativistic particles and the asymptotically free UV behavior can- 
not look like a consistent quantum string theory. Only in the classical limit we can view it as a 
string toy model obeying the classical conformal symmetry. In the classical case it is also easy 
to impose the Virasoro conditions eq.(2.5). In quantum case, we still can try to impose the 
Virasoro conditions, but they will be already some conditions selecting the quantum states. It 
is not easy on the stage of Bethe equations, but we will later make some plausible conjecture 
concerning this selection. 

3.2 The Classical Limit 

In the classical limit, /.i = Cttiq — > 0, the physical mass of the particle mo ~ £~^e^^/^ is 
small in the units of inverse length of the circle on which we put a large number L ^ oo 
of particles, keeping ~ ~ L as well. In this limit the Bethe roots Uk,Vk,6k (and the 
characteristic distances between them) are large, of the order of L. So we can take log^s of 
Bethe equations eq.(3.1) and use the "Coulomb" approximation eq.(3.3) for So{9), as well as for 
other terms: log = + 0{l/u^). ^ Therefore it is useful to rescale 

u = Mx, V = My, e = 

^This limit, sometimes called thermodynamical (which is a bit misleading since the energy is not proportional 
to the length of the system) in the context of ferromagnetic spin chains was suggested in [16] and considerably 
advances in papers [17,18]; see also [24] where an important phenomenon of the anomaly cancelation was observed, 
vital for the validity of this limit. The antisymmetry of the log's of all factors in (3.4-3.6), including log S{9), is 
also important for the validity of the approximation. 
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Figure 1: We plot V{£,) for M = 1,5,9, 13 (lighter to darker gray). It is clear that the potential approaches the 
blue box potential as M ^ oo. 



where we chose 



M 



2tt 



L. 



Taking log of each Bethe equation (3.4-3.6) we arrive at the foUowing system of equations 

1^1 1^1 1^1 



fi sinh ttM^o 



-y— 



+ 



■>■+] ■ 



y 

k 







1 V- 1 2 1 
> + — > 27rni: 



+ 2TTma, (3.10) 
(3.11) 
(3.12) 



where ma,n'j,n1 are integers corresponding to the choice of a branch of logarithm. We will 
call them mode numbers since they are the analogue of the mode numbers of oscillators in the 
string theory. 

It reminds a system of 2D Coulomb charges of different species in static equilibrium, each 
one in its own constant field of force 27rmQ,, 27rn", 27rn^, respectively, and in a confining potential 
for e's 

y(^) = ^cosh(7rM0. (3.13) 

In our limit M ~ log/i~^ — > oo, the potential for ^ ~ 1 basically looks as a box with almost 
vertical walls (the "thickness" of the wall is of the order 1/M) - see figure 1. The positions of 
the walls ^ = it defining the width of distribution of 0's can be estimated with the logarithmic 
accuracy from V{t) ~ 1. One has therefore 



1 

vrM 



(log ^"^ + 0(^1'/°)) ~2. 
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All this tells us that we can write the continuum version of Bethe equations, similarly to the 
large N limit of matrix models, in terms of the resolvents of root distributions: 

1 1 ^ '^^ ^ i i 

13=1 ^'^ 1=1 1=1 ^ 

with the large z asymptotics 

G„(.).Ai, (3.14) 

These equations are defined on the distribution supports Cu^Cq^Cv, each characterized by the 
corresponding integer mode number Uu^riy or m: 

2(f!u{z) - Go{z) = 27rnu, z € Cu 
(pe{z) - G^{z) - Gu[z) = -27rm, z G (3.15) 
2(py{z) — G0{z) = 2T^n^, z e 

There can be as many cuts as there are different mode numbers. Note that we dropped the 
potential term for ^'s in the second eq. since, as we explained, it is zero within the infinite walls 
at ^ = ±2. The only trace left by the potential, as shown in Appendix A, is the behavior of 
Ge{z) ~ ^^j^2 '^^^^ walls of the box which we should impose on the solutions of eq.(3.15). 
Equation (3.15) can also be written in terms of densities defined by 

pe{x)= Yim. Y?y"^(^-^a), 

M^oo IVl ^ — ' 
a 

and similar for u's and u's. Then, for example, second line in eq.(3.15) becomes 

/ P^dw- I P^dw- [ P^dw = -2.m, zeCe. (3.16) 

Jce z-w JCu^-^ Jc. z-w 

We considered only one mode number m. Later it will be shown that in order to reproduce 
the finite gap solutions for classical strings on x Rt we have to assume that we have a single 
cut [—2,2] for ^-distribution characterized with all Bethe roots Ok, k = 1, ■ ■ ■ L having the 
same mode number m. As we will argue in section 5.2 the classical multi-cut states with more 
mode numbers excited correspond to the longitudinal oscillations of the string. In the string 
theory these oscillations are non-physical and forbidden by the Virasoro constraints, but they 
are present in the standard relativistic cj-model eq.(2.18). We know how to do it in the classical 
limit, but it is more difficult to project it in the full quantum space of states. It is also plausible 
to assume that not only in the classical, but also in full quantum theory, only the Bethe states 
with 6a, a = 1, • • • ,L, having the same mode number m obey the Virasoro constraints. We 
don't have for the moment any convincing proof of this conjecture, apart from some examples 
and the comparison of the classical limit of our Bethe equations with the classical finite gap 
solution of section 4. 
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The total momentum can be calculated exactly, before any classical limit 

^ ^ 27r^ sinh(7r6'a) = rriL - ^ UpSp - ^ UpSp (3.17) 

a p p 

where Sp, Sp are the filling fractions, or the numbers of Bethe roots with a given mode numbers 
n-u,p,ny^p. To prove it it suffices to take the sum of logarithms of the eq.(3.4) for all roots 6^- 
The contribution of Sq{9) terms cancels due to antisymmetry while the second and third sums 
in the r.h.s. of (3.10) are excluded using (3.5) and (3.6) respectively. 

For the closed string theory we should take P = 0, which gives the level matching condition. 
For the perturbative super YM applications one should take 5^ = [69]. Then we have the well 
known formula "^pnpSp = uiL (see [19] for details). 

Let us now study the system of equations eq.(3.15) as a Riemann-Hilbert problem, first for 
the highest weight states with U (1) symmetry, without magnon excitations Uk and v^, and than 
for a general classical state. Our goal will be to reproduce (and to generalize from a single sector 
to both left and right sectors) the finite gap results and the corresponding algebraic curve of 
the paper [19]. 

3.3 Highest Weight States of U{1) Sector in Classical Limit 

If the right and left modes are not excited we have only the states with U (1) modes, which 
means that the currents ji = are diagonal matrices. 

In the classical limit, using the Coulomb approximation, we have for this sector the following 
Bethe equation 

1 ^ 1 
u sinh vrMf Q, — 27r m = — — > t 

In the continuous limit, the equation for the asymptotic density, L ~ M — > oo, is given by 

^eix) = -2-nm, x e Co , (3.18) 

with inverse square root boundary conditions near ±2 (see Appendix A and discussion in previ- 
ous subsection). The asymptotic density of rapidities, pe{z), is obtained through the resolvent 
from 

Pe{z) = -77^ {Ge{z + iO) - Ge{z - iO)) . 
2m 

The analytical function Gq{x) having a real part on the cut defined by eq.(3.18), with support 
[—2,2], with inverse square root boundary conditions and behaving at z — > oo as ^^(-z) — > 
is completely fixed: 

G0{z)= i—^==JL-27rm\ (3.19) 
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Figure 2: Density of S-roots before and after phase transition. Black line - asymptotic densities of eqs.(3.20,3.28), 
blue dots - numerical solution for L — 150 roots. 

which gives for the density 

if 27rm z + jj\ , ^ 

Notice that the distribution has a singular behavior at the endpoints which will be the typical 
behavior even for the general multi-cut solution which we study in the next subsection. Notice 
also that applying to the eq.(3.19) the Zhukovsky map 

1 

z = x + - 3.21 

X 

we obtain 

Geizix)) = ^ — + , ^ (3.22) 

X — 1 X + 1 

which shows the poles at x = ±1, typical for the finite gap solution of the section 1. The 
Zhukovsky map will be the central piece of our proof of the identification of the continuous limit 
of Bethe ansatz equations with general classical solutions of cr-models considered in this paper. 
From the general formula eq.(3.17) the momentum of such a state is 

P = ml (3.23) 

Let us now compute the energy of this state in the continuous limit (it is derived in greater 
detail in Section 3.5 for a general state). We have to compute the sum 

E=—y^ cosh(7rM^„) ~ V ^ sinh (vr^^M) - V ^ sinh (nCaM) 

" ^a>0 5a <0 

where the last equality holds with exponential precision since the decaying exponents of the 
sinh^es can be dropped. Then, by the use of Bethe ansatz equations for each of the sink's we 
get 

- V log Si (M [e„ - + V m sign(C„) . 
vr ^ — ' ^ — ' 

C/3<o<Cq 
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Now we do not have a problem, as in eq.(3.24), with badly defined sums of small exponentials. 
In the scaling limit these sums become integrals 



E ^ -— [ dz [ dw^^^^^^^ + mM [ peiz)sign{z)dz (3.24) 

TT J-2 Jo Z-W J 

where peiz) is given by eq.(3.20). In fact these integrals can be easily calculated and the result 
reads 

7-2 

E = —— + 47rMm2 . (3.25) 

This calculation is generalized to an arbitrary solution of eq.(3.4)-(3.6) including the left and 
right excitations (with arbitrary u^s and v^s) in Section 3.5. 

Now notice that if we start decreasing the number of particles L, at the point 

L* = A7T\m\M (3.26) 

the density (3.20) will vanish at one of the walls. This point certainly signifies a phase transition 
in the effective spin chain and a non-analytical change in the classical state. If we continue to 
decrease L past this point the new support will now stretch only on the interval [—a, 2] (for 
m > 0) and the distribution in the limit z — > —a becomes less singular: pe{z) ~ \/a + z. The 
resolvent satisfying the right behavior at infinity is: 



2+ ^ 



Ge{z) = 2vrm | ' ^j^^^' " ^ I (3-27) 



which gives for the density 



p,(z)=2m^/|±| (3.28) 



where a = — 2 + r^f^jj- The phase transition corresponds to a = 2. This phase transition fits 
nicely with the electrostatic picture point of view since m plays the role of an electric field 
which pushes the particles against one of the walls. It is also interesting to note that and 
K_ become of the opposite signs, which, according to 2.4, means that the coordinate Y becomes 
space like on the world sheet. 

Only L > L* looks compatible with the description of the classical finite gap solution of 
section 2 for the SU{2)l x SU{2)r sigma model. Only in that case the poles in x plane will 
be on the right place at x = ±1. For the string theory the total momentum should be zero 
P = 27rm = 0, hence the classical limit works well, at least in this highest weight sector. In the 
presence of L,R-excitations with non-trivial Ui^Vi roots the phase transition is still possible. It 
would be very interesting to understand what are the states beyond this phase transition and 
what classical solutions correspond to such states if any. 
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3.4 General Classical States and Its Algebraic Curve 

Let us now reproduce the general classical states in terms of an algebraic curve as a solution of 
the Bethe ansatz equations (3.15) in the scaling limit. 
We define the quasi-momenta by the formulas 

Pi = -P2 = Gu- ]^Gg 

P3 = -P4 = G,-^Gg. (3.29) 

Then the functions p'i{z),P2{z) have the cuts of the type C^, Cg, and p'^{z) , p'^{z) have the cuts 
of the type C^^Cg. As in the previous section, we consider only the situation with a single cut 
Cg and any number of C^, Cy cuts. We notice that the quasi-momenta p'^{z),p'2{z),p'.^{z)^p'^{z) 
form four sheets of the Riemann surface of an analytical function p'{z) (see fig. 3). Indeed 
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Integrating these equations we restore the mode numbers n„, and m, 



X G Cm, 


PI-" 


-P2~ 


= 27rn„ 


X G Ce, 




-P3^ 


= 27rm 






-P4~ 


= 27rn„ 


X G Cg, 




-pr 


= 27rm 



(3.30) 



These equations can also be written as holomorphic integrals around the infinite B-cycles: 

dp = 2TTnu,j rij = 1, . . . ,Ku 

dp = 2-Knyj rij = 1, . . . , (3.31) 
dp = 27rm 



where the the first two conditions correspond to the equations in the first and third line of 
(3.30), respectively, while the last one corresponds to any of the equations of the second and 
fourth lines of (3.30). The B cycles are defined as in fig. 3. 

From (3.14,3.7) we have the following large z — > oo asymptotics of quasi-momenta on differ- 
ent sheets 

Psi^)^-^h P^(^)^JMh (3-32) 
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The fiUing fractions, or the numbers of Bethe u, f-roots forming each cut, are defined as 
follows: 

M r M r 

S'l = ir^i> P{z)dz, Sr = —f p{z)dz (3.33) 

These contours are represented in figure 3. 

3.5 Energy and Momentum 

In this section we will show that the constants k± defined as the square root singularities of the 
quasi-momenta 

pi,3(z) = T^==, \A > 2 (3.34) 
\/±z -2 

are completely fixed by the energy and momentum. At this point we know that (3.17) 

P=^Y1 smh(7r^„) =mL-Y, npS; - ^ nj,S; . (3.35) 

a p p 

Let us now calculate the energy E. As a byproduct we will also compute the momentum in 
terms of the singularities at z = ±2 described by k,±. 
We want to compute the sum 

^ = ^ J^cosh(7r0„), 

a 

but we cannot simply replace this sum by an integral and use the asymptotic density for ^'s 
to compute the energy. This is because the main contribution for the energy comes from large 
^'s, near the walls, where the expression for the asymptotic density is no longer accurate. It is 
natural for the classical limit since the particles become effectively massless, the contributions 
of right and left modes is clearly distinguishable and are located far from 9 = 0. 

The calculation (3.35) of momentum was simple because we had to sum over sinhvr^a which 
can be taken from the Bethe equations 

/isinhvr^^ = iY^logS^{da-e/3)-Y^-^ a \ + 27r?n . (3.36) 

Now we have to sum over coshTr^Q. We notice that the energy is dominated by large 0's 
where, with exponential precision, we can replace coshvr^Q, by ztsinh-Tr^Q for positive (negative) 
9a- Furthermore the contribution from the 0's in the middle of the box is also exponentially 
suppressed since fi is very small. Thus we can pick a point a somewhere in the box not too close 
to the walls. One can think of a as being somewhere in the middle. Then, 

E= - sinh (TT^aM) - y ^ sinh (vr^aM) , 
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where, let us stress, the result is correct and independent of the point a with exponential precision. 
Having a sum of sinhvr^a we can substitute each of them by the corresponding Bethe equation 
(3.36) obtaining 

7-, \ ^ , ^9/,^r.- A IN \ ^ sign (fa — a) \ - signffa — a) \ - . ,^ , 

" =^ i<?<r '° - - £ - ^ - ? - - ■ 

Now we can safely go to the continuous limit since in the first term the distances between ,^'s 
are now mostly of order the 1. Moreover, it is very important that the contribution from the 
^'s near the walls ±2 is now suppressed since 

I log Si{M{2 -^f3)\>\ log Si{M{2 - a)| ~ 1/M . 

This allows to rewrite the energy, with 1/M precision, as follows 

E ~ / dz \ dw / signfz - a) dz dw 

■K J_2 Ja Z-W 2tT J Z-W 



M f P0(z)py(w) . , s , , r f / \ ■ r \ > 

— / sign(z — a) dz dw + mM / p0(zjsign(z — a) dz 

2tt J z — w J 



(3.37) 



where the densities are the asymptotic densities for the exact box potential, given by the integral 
eq.(3.16). 

As a result, by the use of Bethe equations, we managed to transform the original sum over 
cosh, highly peaked at the walls, into a much smoother sum where the main contribution is now 
softly distributed along the bulk and where the continuous limit does not look suspicious. 

Let us now take (3.37) as our starting point. From the previous discussion we know that 
this expression does not depend on a provided a is not too close to the walls. In fact, we can 
easily see that it does not depend on a at all after taking the continuous limit leading to the 
perfect box-like potential. To see this one can see that, due to Bethe equations eq.(3.16), the 
a-derivative of this expression is zero for all a g] — 2, 2[. Hence we can even send a close to the 
wall: a = —2 + e, where e is very small. Let us calculate the first term. The main contribution 
to the integral comes from — 2 ~ z ~ if so that we can use the asymptotics (3.34) to get 

.-2+e .2 pg{z)pe{w) f~^+' , , 4Mk2 

/ dz / dw ^^ ' ~ - / dz \ dw , ~ 2ttMk_ 

Vr J~2+e Z-W J ~2+e T^{z - w)\J 2 + Z^j2 ^ W 

The remaining 3 terms are very simple: since a ~ —2 we can simply drop the sif/n-functions 
inside the integrals and obtain exactly the expression of the momentum in the continuous limit. 
We arrive therefore at 

E ~ 2Mk^^ + [ mL - ^ nlSl - ^ n^S^ j . (3.38) 
\ V V ) 

where the expression in the parentheses is the momentum. If we compute the a-independent 
integral (3.37) near the other wall, i.e. for a = 2 — e, we find 



E ~ 2M 



\ V p J 
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Figure 3: Structure of the curve coming from the Bethe ansatz side. The quasi momenta pi,2,3,i{z) are defined 
in (3.29). This figure is related with fig.4 by means of Zhukovsky map. 

Therefore, equating the results one obtains the desired expressions for the energy and momentum 

E±P = 2ttMk^^ (3.39) 
through the data k± at the singularities of the curve at z = ±2. 

4 Matching with the Finite Gap KMMZ Solution 

In this section we will demonstrate the main result of our paper: the equivalence of the large 
density limit of a system of physical particles on a circle for the quantum 0(4) cr-model described 
in the previous section, to the classical theory of the same model described in section 2. On the 
one hand we have the quantization of the sigma model given by a set of Bethe equations. For this 
quantum system Virasoro constraints are somewhat obscure. We introduce however a natural 
selection rule for the quantum states, such that in the scaling limit we obtain the classical finite 
gap KMMZ solutions for the x classical sigma model with classical Virasoro constraints 
imposed [19]. Let us also recall that this sector belongs to the superstring theory on AdS^ x 
and is in itself a consistent classical truncation. 

Our result also means that the considered scaling limit of the quantum 0(4) cr-model is 
nothing but the classical limit and our result gives a decisive demonstration of the correctness 
of Zamolodchikovs' physical bootstrap S-matrix approach [38] to these models, as well as of the 
Polyakov-Wiegmann solution of the principle chiral field [39]^. According to our results, both 
approaches construct indeed the quantization procedure with the correct classical limit. ^ 

More precisely, we will show here that every solution of Bethe ansatz equations in the 
scaling limit, described by an algebraic curve of the quasi-momentum (3.4-3.6), is in fact a 

**Altliough we demonstrate it here only for the SU{2) case we see no obstacles for the generalization of our 
method to any SU (N) principle chiral field. 

''Modulo the phase transition noticed in the previous section, to the states which are hard to identify with the 
classical solutions of the underlying classical model. 
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Figure 4: Algebraic curve from the Finite gap method of KMMZ. In this language, u and v cuts correspond to 
cuts inside and outside the unit circle respectively. This figure is related with fig. 3 by means of Zhukovsky map. 



finite gap solution for tlie classical string on x Rt. In the next section this correspondence 
will be generalized to the 0(6) sigma model corresponding to the whole x Rt sector of the 
classical superstring. The generalization to all 0{2n) sigma models is also straightforward (the 
corresponding classical string, unlike the quantum one, is well defined). 

To prove this correspondence, we have just to compare the algebraic curves describing a 
Bethe state in the scaling limit and the finite gap curve eq.(2.30), together with all their moduli 
and the data at singularities. The similar goal for the OSP{m + 2n\m) model was achieved 
in [64] by the direct solution of integral Bethe equations. 

Central to this comparison will be the Zhukovsky map z = x + 1/x. We will show that 
the Riemann surface in z variable of the scaled Bethe equations, fig. 3, on the one hand, and 
the Riemann surface in x variable of the classical finite gap solution eq.(2.30), fig. 4, on the 
other hand, are two different projections of the same algebraic curve related by Zhukovsky map. 
We attempted to present schematically two different projections by different colors of projected 
parts of the curve fig. 5. 

When we apply Zhukovsky map to the Riemann surface of fig. 3 the 0-cuts along (—2,2) 
disappear on all 4 sheets and these singular branch points become simple poles, since under the 
map 

, ^ < — • (4.1) 

^/I±2 X±l ^ ' 

The Riemann surface in the x projection will consist only of two sheets, as on fig. 4. All the 
u-cuts and u-cuts connect now these two sheets. 
The inverse map 

x± = ]^[z± ^/z^ - 4) (4.2) 

projects the two upper sheets of fig. 3 with n-cuts into the interior of the unite circle in x 
projection, by means of x_(z), where as the two lower cuts of fig. 3 with t;-sheets are projected 
into the exterior of this unit circle on x projection, by means of Xj^{z). More precisely, we will 
show that 

P^{z) = p{x+{z)), pi{z) = 2TTm - p{x-{z)) . (4.3) 
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Since, by definition, p2 = —pi and p4 = —p^ we see that ^3^4 describe the exterior of the unit 
circle in the upper (lower) sheet while p2,i describe the interior of the unit circle in the upper 
(lower) sheet. 

Comparing the singularities at x = ±1 after the map we see that they are again given, due 
to eq.(3.34), by the energy and momentum in terms of k± in full correspondence with the finite 
gap result eq.(2.28). 

Under this map the 4 infinities of the 4 sheets of quasi-momenta pi,P2,P3,P4:, fig-3, are 
mapped according to (see (4.2,4.3)) 

(002,1,003,4) i — > (0±,oo±) . 

Then, from the first line of (3.32), we read off the behavior of the quasi-momenta near x = 0± 
while from the second line of (3.32) we infer the behavior of the quasi-momenta at x = oo-t. 
They match the asymptotics (2.27) and (2.26) respectively after imposing M = ^ or 

-log/U = ^, (4.4) 

a simple and nice relation exhibiting the dimensional transmutation property of the asymptot- 
ically free theory with the precise coefficient^'^. 

The filling fractions (3.33) coincide with (2.34) of the finite gap solution since they are 
homomorphic integrals, invariant w.r.t. the change of projection^^ . The same can be said about 
the Bethe equations eq.(3.31)) which we have to compare with the finite gap equations (2.33, 
with the following definition of the mode numbers 

= 2m-ni , n\ = ni. (4.5) 

^"For the 0{N) sigma model the beta function for the couphng is given by /3 = g-f^^ \/ •^(^) = N — 2 where A 
is the cutoff of the theory. The dynamically generated mass must be of the form m = A/(\/A). The functional 
form of / is fixed by the /3 function upon imposing independence on the cutoff of this physical quantity. Thus, 
for general iV, -log^ = ^ + 0(1). 

^^The minus sign in (2.34) comes from relation (4.3) between pi{z) and p{x). 
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This redefinition is due to the fact that a -B^-cycles of eq.(3.32) going from ooi on the upper sheet 
of fig. 3 through the u-cut to 002 on the next sheet, are projected by Zhukovsky map into the 
interior of the unit circle of the x projection, fig. 4, with the 001^2 projected to 0^=. To complete 
it to the S-cycles as defined in eq.(2.33) we have to complete them with paths (oo_,0_) and 
(0+,oo_|_), each of them is equal to quasi-momentum —p{0) = — 27rm. This completed cycle will 
run in the opposite direction of a usual cycle, hence the minus sign in eq.(4.5). 

The third equation eq.(3.31) follows immediately since integral over the iJg-cycle is mapped 
into the integral from cx)+ to 0+ yielding p{0) as result. 

We conclude that the two projections, one obtained from the quantum Bethe ansatz, the 
other from the classical finite gap solution, represent the same algebraic curve with the same 
moduli and thus the underlying solutions (states) completely coincide. 

As a result we can express densities of the Bethe roots thought classical quasi-momentum 
p{x) as^^ 



Pv{z{x)) 
Puiz{x)) 

g{2 cos (p) 



p{x), xeC" 

pix), X G 

_lim f^^'^ 

vr 



G [0,7r] 



where 2TTi p{x) = p{x — iO) — p{x + iO). These densities automatically satisfy Bethe ansatz 
equations eq.(3.16). 

In Appendix B we provide an alternative, completely algebraic, derivation of the equivalence 
between the curves. 



5 Some Limiting Solutions of the Quantum 0(n) Sigma-Model 

Here we will consider three limiting cases: the BMN solutions from the quantum Bethe equa- 
tions, the semi-classical solutions with small longitudinal amplitudes and the discrete spin chain 
limit. 



5.1 BMN Limit 

In this section we will consider the BMN limit [72] which, in our notations, reads M ~ L ^ 

1, J^, Jm ~ 1. We are therefore considering a big 9 cut which will be slightly perturbed by some 

point like microscopic u and v cuts, treated as corrections. The position of these microscopic 

cuts in the presence of the 9 cut are found from the first and third equations in (3.15) with Ge 

given by (3.19) with m = (the only choice compatible with the stringy condition P = for 

^■^Last expression deserves some words. Points 2 ± iO below and above the 6 cut are mapped to a; = x+{z) and 
X = X-{z) = l/x+{z) Since z G [—2,2] these are conjugate points in the unit circle. Then the discontinuity of 
p{z) in the 9 cut is given by the imaginary part of p{x). 
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very small filling fractions for n's and -y's.). We find 



LI / L2 

27rn, Xn = sie:nn\ A + 



Then, from the second equation in (3.15), we can compute the correction 5Gg to the resolvent 
of the main cut. 



^ X Xn 

where Nn = Jn/M. Since 

6Gb{x) = 5(p0{x) + i7T6p0{x) 

we have 

X r \ ST AT ( \/a;2 -4 signx^ 1 

^ \ 7r\/4 - x2 x-Xn 7rV4 - x2 



(5.6) 



Let us explain why. The inverse square roots follows from the known support and behavior of 
p near the walls. Then the first term inside the brackets follows from cancellation of the poles 
at Xn in SGq. Finally the second term in fixed by requiring 0(l/a;2) decay of 6Gg at infinity, 
or, equivalently, from requiring that the perturbation of the density does not change number of 
particles in the cut, / 5pg = 0. 

From (5.6) we can read the change in the behavior near ±2 of the density of rapidities, i.e. 
compute 6k± and, therefore, 

5P = 27rMK{5K+-5K^) = -^nJ„ (5.7) 

6E = 2t:Mk{5k+ +5k_ ) = A ^ j„ _ ^1 + . (5.8) 

The same result was obtained in [19], but in a different regime, the semi-classical one, where 
K J„, J„ < L ~ M. In our case are small intcgcvs which refiects the quantum nature of 
the result, as for the original BMN result [72]. 

5.2 Multi Cut Vacuum States 

In this section we will try to understand the meaning of the multi-cut solutions in 9 and its 
possible relation with Virasoro constraint. 

We will consider here the solutions in the Abelian U (1) sector described only by ^-variables: 



9 = exp 



-ITS 



irr + Sa + ^A^ cos(w(r + <r)) + Y^ A^ cos{uj{t - a)) j 

\ u>>0 uXO / 



(5.9) 



where S is the winding number, integer to ensure the periodicity in a, and T = Qr/VX is 
the angular momentum. For small amplitudes A^ we can quantize this solution as a set of 
oscillators. The energy in this approximation reads 
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where N^^ is the mode number of the oscihator with the frequency uj. This formula is a well 
known expression for the spectrum of U (1) model. It is very natural since we can use arguments 
of BMN to see that in the limit of large angular momentum the system feels only a small 
neighborhood of the main circle on and the longitudinal oscillations are described by a 
massless scalar field. This is why our analysis in only valid for small amplitudes Am- 
The worldsheet momentum is 

P = VXST + 0{1). (5.11) 

The terms 0(1) represent quantum corrections which are not captured by the leading semi- 
classical approximation. 

We will show that the solution eq.(5.9) with small amplitudes corresponds to a multi-cut 
configuration of 0's when all cuts except one are small (with filling fractions of the order 1). We 
will denote by tuq the mode number of the large cut and by L^o the number of 0's in this cut. 
We assume Lq ~ \/A. 

We can compute the momentum and the charge in terms of the Bethe ansatz quantities 

P = nioLmo + ^ mLm, Qr = L^o + ^ Lm = L (5.12) 

and thus 

r=A, s = mo (5.13) 

where we are allowed to drop 1/Lq I/'v/A terms due to the winding number quantization, 
5 G Z. Having identified S and T in the language of Bethe ansatz we can write 

where mo and L are integers. It is left to identify the amplitudes with the remaining filling 
fraction. Let us denote A''^^ = Lcj+mo 1 then we have 

£' = £'0-1- ^ \m-mQ\Lm- (5.15) 

m^mo 

Expression (5.14) can also be established form the Bethe ansatz side. One has, see section (3.3), 
equation (3.24) 



E = - log Si {M[io.-ip\) + Y.m^^\gii{i^). (5.16) 

C/3<0<Cq a 



The main cut with the mode number mo occupies most of the box. Other cuts are squeezed 
close to the right (left) wall for m^ greater (smaller) than mo- Then the second term in (5.16) 
reads 

^ rria sign(^a) = X] "^0 sign(^c.) + Y\m- mo| (5.17) 
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We only have to show that the change in the first term in (5.16) and in the first term in (5.17) 
is smah compared to the last term of (5.17). We wih demonstrate that the change in one of the 
mode numbers leads to 1/L displacements of ^'s. Consider for example the situation when we 
decrease the mode number of the first by (5n ~ 1. This 6 has the biggest displacement since 
the force acting on it increases by lirdn. All other ^'s are affected only due to the displacement 
of the first one. Let us estimate this displacement assuming it to be small. Then one has the 
following equilibrium condition 

7r/icosh(7r6li)(56'i ~ 1. (5.18) 

Since vr^u cosh(7r^i) S> 1^'^ the displacement is small 59 <^ 1, then 6^ <^ 1/L. This means that 
the change of the density is negligible 5p(C) <^l/L and we can say that the first term in (5.16) 
remains intact. 

We obtain the precise agreement between semi-classics and Bethe ansatz calculations for 
several 9 cuts. We can conclude that the longitudinal oscillations, which break the Virasoro 
conditions and thus being inadmissible in the string context, manifest themselves as extra cuts 
in 9. We observe indeed this phenomenon, at least when these oscillations are small. As it 
was said already in the introduction, we conjecture on the base of these arguments and of the 
comparison of our classical limit with the direct finite gap approach that the Virasoro conditions 
impose the selection rule on the states described by (3.4,3.5) under which only the solutions 
having ^^'s with the same mode numbers mi = ■ ■ ■ = mi are retained. 



5.3 XXX Spin Chain Limit 

In this section we will consider the XXX spin chain limit of the Bethe equations for the principal 
chiral field. Contrary to the classical limit considered above it is rather a strong coupling limit. 
Namely, we take ;U — > cxo and arbitrary L, J„, J^. In this limit ^'s are squeezed near zero. In the 
leading order we can just take them equal to zero so that (3.5,3.6) become non-interacting two 
SU{2) XXX chains 

= ]T "^'""' + \ (5.19) 

Uj — i/2 I Ui — Uj — i 

Vk + i/2\^ YjVk-vi + i 



n 'I _: • (5.20) 



Vk-i/2j f-_^]:'"k-vi 

These equations look similar to the = 4 SYM XXX spin chain Bethe equations, the first 
one for the sector of X, Y scalars and the second for X, Y scalars. To make a more precise 
comparison we have to fix the conserved quantities. Let us firstintroduce the notations 

u + is/2 



es{u) = i log 



u — is/2 



^^We know that E = j^/iX^a coshTr^a ~ L. With exponential precision we can sum only up Oa < —2AI + L^^^ 
so that number of terms is of order L^'^* as one easily sees from the density computed in 3.20 . Then we have 
E < fi cosh{n9i)L^^'^ which means ^cosh(7r^i) ^ 1 
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and 



Below we will show that energy of the system given by eq.(3.9) is proportional to that of the 
XXX spin chain 

To see this we expand further in powers of /u. The formula for energy with 0{9'^) precision 
becomes 



-£eG.^). 

Q = l ^ ^ 



(5.22) 



Let us evaluate 9^ and 9'^. We expand BAE (3.4) in powers of 9 up to 0{9'^) 

L 

/xsinhvr^Q, ~ 27r m + (tt sign(6'Q, — g/j) — A{9a — 9p) log 2) 



(5.23) 



k j 

The second term in the r.h.s comes from the expansion of Sq. Expanding eqs. (3.5,3.6) we have 
27m^ = Lei(nj) -L ^e2('Uj- --Ui) -eU%) J^^a + 0(^^), (5.24) 

27rnl = Lei{vk)-LYe2ivk-vi)-e[{vk)Y^» + 0{9^). (5.25) 

l^k a 

Summing all eqs. (5.24,5.25), using the level matching condition P = and eq.(5.23) one obtains 



/ivr^Q, ~ 7r(2a — L — 1) + 



9/3 {E^..-4Llog2) + 0{9^), 



so that 



90 = 2a - L - 1, 



0i = 9'i (Sx^-4Llog2) 



and thus, from eq.(5.22), 

"L/i 7rL(L2 - 1) 27rL2(L2 - 1) 



E 



2tt 



+ 



log 2 



+ 



^L(L2 - 1) 



(5.26) 



(5.27) 



12/x 3/i2 

This formula, together with eq.(5.19), reproduces the energy and the Bethe equations for the 
XXX spin chain up to the constant terms inside the square brackets. Although it is known 
that the perturbative expansion of the energy E in general does not coincide with that of the 
eigenvalue of the dilatation operator in SYM [28] , it has been confirmed that they coincide up 
to two-loop approximation in the continuous limit [19]. The above result may serve as a starting 
point to elucidate the matching at the discrete level. 
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6 0(6) Sigma-Model 

Let us now move on to a larger subsector of the superstring theory, namely let us consider the 
0(6) non- linear sigma model. As before we describe a general quantum state by the system of 
vector particles which live on a circle of length C = 2tt. The wave function is then specified by 
a set of rapidities {Oa} which determine the coordinate part of the wave function, plus a set of 
Bethe roots {u^^\u^^\u^p} encoding its color structure^^. When imposing periodicity of the 
wave function one obtains the quantization of these parameters. They are are constrained by 
Bethe ansatz equations (see Appendix D for the derivation and for the generalization to 0{2n)) 

^ ft 7/^^) -^779 

1 = n 

L (2) n K2 „,(2) (2) (3) „■ /n Xi (2) (1) . 

7/(2) ft 772 ^^iP^ iP^ i^^iP^ ,/3) +,72 -'■^7/2) _,,(!) ,„72 ^ ' 

1 = n 

where 

r(i-.f)r(i-zf)r(| + ^f)r(i + ^g) 
°^^""r(| + .|)r(i + .f)r(|-^f)r(i-.f)- ^'''^ 

In the next subsection we will show that, in the classical scaling limit L ~ log f ^ 00, the 
solutions of this system of Bethe ansatz equations are in one to one correspondence with classical 
solutions of 0(6) sigma model, classified by means of an algebraic curve [21]. The derivation 
here will be similar to the SU{2) principal chiral field considered in previous section. In the 
classical limit when ft's will be large we can use Coulomb approximation and substitute Sq{6) 
by its large ft asymptotics 

Hog So{ft) = ^ + 0(1/0^). 

Furthermore, the potential becomes a square box and Oa/M will therefore be distributed, at 
leading order, between —2 and 2 provided we define the variable z = ft/M with M = — 
Thus, in this limit, we can recast (6.1) as 

2(^1 -G2 = 27rn(^) 
2(^2 - Gi - G3 - Ge = 27rn(2), (^g - G2 = 2TTm (6.3) 

2(^3 -G2 = 27rn(^) 



^''For a general simple group one would have r different kinds of roots where r is the rank of the group, i.e. the 
number of simple roots. 
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where the resolvents are given by 

1 1 I ^'^ 1 
Ge(z) = — V — — , Gi(z) = —y , / = !,..., 3, (6.4) 

M^z-Oo^/M' ^ M^z-u(')/M ' ' ' \ J 

and each equation holds for z belonging to a or n^*^ cut if Ge or Gi is slashed, respectively. 
These equations, as we will show in the next subsection, can be used to define some algebraic 
curve which maps onto the classical curve of [21] after Z-transformation z = x + 1/x. 

The quantum state (Bethe vector) corresponding to the solutions of (6.1) carries the following 
so (6) spins 

Jl = L-K2, 

J2 = K2-Ki-K^, (6.5) 
J3 = Ki-K^, 

where (Ji, J2, J3) are measured by the orthonormal basis of the Cartan subalgebra. 
6.1 Algebraic Curve 

To relate solutions of BAE (6.1) with classical solutions classified by algebraic curves [21] we 
will construct possible curves corresponding to (6.3). This means that the resolvents Gi^2,2,{z) 
and Ge{z), when taken in appropriate linear combinations, are in fact different branches of 
some unique analytical function. In fact there are several possible combinations, each of them 
corresponding to different representation of 50(6). 

6.1.1 Curve for the Vector Representation 

The curve in vector representation is the most simple one. Introducing quasi-momenta 

qi = Ge — G2 

q2 = G2-G3-G1 (6.6) 
Q3 = Gi — G3 

one can easily see from eq.(6.3) that ((/i 2 3' ~1i 2 3) branches of the same function. Riemann 
surface corresponding to this function is depicted in fig. 6, to the right. In terms of these quasi- 
momenta eq.(6.3) takes the form 

^1-^2 = 2^2) 

^s-h = ^^ri^^\ ^i = 27rm (6.7) 

h+k = -27rn(3) 
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Figure 6: Riemann surfaces resulting from classics and from the Betlie ansatz match after Zhukovsky map. On 
the left picture we plot some possible cuts coming from the classical curves analysis, (2.17). For each cut outside 
the unit circle there is a mirror cut inside the unit circle, (6.10). In this picture we did not plot some possible 
cuts, for sake of simplicity. There is no cut, for instance connecting gi and 53. On the right picture we present 
the curve obtained from the Bethe ansatz point of view, (6.7). Here it might seem that we have plotted every 
possible type of cuts. It is not the case - due to the existence of stacks [24] roots of different kinds can attract 
each other forming some extra cuts. This extra cuts would be the image, under Zhukovsky map, of the cuts which 
we did not plot on the left figure. 



Let US see how this curve appears from the results of [21], reviewed in Section 2.2. There 
one has the curve iqi,q2,Q3, —Qi^ —Q2^ "^s) with the properties (2.14,2.15,2.16), 

q,{x) = 6k=i^^ + 0{l), (6.8) 

Ux) = ^- + 0(l/x2), (6.9) 

qk{l/x) = (47r?n - gi(x)) + (1 - (5fc,i)gfe(x) . (6.10) 

Charges Jj are properly normalized to be integers. 

Let us check that Z-transformation maps this classical curve to the curve (6.6) corresponding 
to a solution of BAE (6.1). Namely^^ 



qk{z) = qk{x+) x+ = -yz+ - . (6.11) 

Where the square root is defined in such a way that |x+| > 1 for z G C. Since the values of 
qi inside the unit circle are related with the values of qi outside unit circle by eq.(6.10) we can 
always reconstruct q{x) in the whole x plane from q{z). Prom (6.8,6.10), one has, for small e, 

00 

(?i(l + 6) = 2vrm + ^ + X]ci'^e''-' (6.12) 



fc=i 



^2,3(1 + e) = E^f'^^'' (6-13) 

fc=0 



^^We notice that this is the only possible normalization ~ (6.7) and (2.17) forbid us to multiply qu by any 
constant coefHcient. 



32 



Since, after Zhukovsky map, qk{2 + e) = ^^(1 + y^), one sees that a new cut appears connecting 
^q'l while no singularity appears for g'j 3 (see fig. 6). This is consistent with eq.(6.6) where only 
qi has 9 cut. 

Let us now understand this Z-map in greater detail. We start with six sheets (±g^ 2 3)- Then 
we apply Zhukovsky map by means of which each of these sheets is mapped to a new pair of 
sheets, one coming from points located originally inside the unit circle while the other is the 
map of the points in the exterior. We designate this set of sheets by {^qi2z^^l'i2\)- Now, 
from (6.10), we know that, 

(Tgr*,T92r3"*) = (±9'r,T92:3)- 

Thus we can keep only the sheets in the left hand side. Their cut structure, figure 6, is inherited 
from the original ones with an additional cut between itg^ as explained above. It is exactly the 
structure we find from the BAE point of view (6.6). 

A small remark is in order here. From (2.17) we see that classical curve can have more cuts 
than are listed in eq.(6.7) and depicted in fig. 6. For example and q'^ can be connected by a cut. 
This apparent discrepancy is solved by the introduction of stacks [24] . In the thermodynamical 
limit Bethe roots of different types can group near the same extremum and form cuts. To 
describe them one can also define resolvents. As a result, quasi-momenta eq.(6.6) can have cuts 
connecting for example qi and (73. For more details see [24], Section 4. 

It is now clear, that the equations in the first row in (6.7) follow from (2.17) using the map 
(6.11). It is also the case for the equation to the right row. Indeed, from eq.(6.10), one has 

2^i{z) = qi{z + m) + qi{z - m) = qi{x+[z + iQ]) + qi{x+[z - m]) = A-nm (6.14) 

where for the last equality we use Xj^{z + iO) = 1/2;+ (z — iO). Other equalities are just restate- 
ments of our definitions. 

From (6.12,6.13) and (6.6) we see that 

Pe = -^{Ge{z + i{))-Ge{z-iQ))c^^=^, z ^ ±2 (6.15) 
2m \J2^ z 

Now, having related q\^2,z{.z) and qx^i^ix) we can relate \f\ and M. Namely from the 
definitions (6.6,6.4) we immediately see that, for large z, 

\L-K2 



91(2) 



z M 



1 K2-K^- K\ 
- , (6.16) 



IKi-K; 



3 



z M 

Comparing with (6.9,6.5) we obtain M = or, since M = — 



logM = ^. (6.17) 
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Figure 7: Classical curve from [21]. For each cut outside the unit circle there is a mirror cut inside the unit circle 
- eq.(6.22). In each sheet there are poles at ±1 which will lead to 6 cuts after Z-transformation - see fig. 8. 



This result matches perfectly the 1-loop /3-function of 0(6) sigma model (see footnote in Section 
4). It is a strong confirmation of our correspondence between the quantum Bethe ansatz eq.(6.1) 
in the scaling limit and the classical sigma model. 

6.1.2 Curve for the Spinor Representation 

In this subsection we will construct the curve associated with the spinor representation. Having 
already seen the mechanism at work in the previous section we will now proceed rather briefly. 
An 8-sheet surface can be build by introducing quasi-momenta 

Pi = ^Gg-Gi 

P2 = \Ge + Gi-G2 (6.18) 
P3 = -\Gg-Gz + G2 

PA = -^Gg + Gs. 

It is then easy to see that (p'^ 2 3 4' ~P'i 234) branches of the same analytical function (see 
fig. 8). Indeed 



Pi -Ale, = 


-2(^[ + G2 = 


Pi+P2\ce = 


(^'g -G[=0 


P2-P3\c, = 


-2(^'2 +G'g+ G[ + G'3 = 


p'3+p'a\co = 




P3-p'a\c, = 


-2(;'3 + G'2 = 



(6.19) 



This curve corresponds to the classical curve in spinor representation [21]. As in the previous 
subsection they are related by Z-transformation (see fig. 7). After Zhukovsky transformation 4 
sheets split into 8. Since there are poles at ±1 in each of the original 4 sheets the resulting 8 
sheets will be pairwise connected by 4 cuts as shown in fig. 8. 
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Figure 8: Curve corresponding to the classical limit of the quantum Bethe ansatz. It is related by Z-transformation 
with the curve on fig. 7. 

Let us use the spinor curves to relate M and \/A as in previous subsection. The vector 
representation (of the previous subsection) and the spinor representation are related by [21] 

gi + 92 - 93 9i - 92 + 93 
PI = ^ , P2 = ^ 

-9i + 92 + 93 - 9i - 92 - 93 

P3 = ^ , P4 = ^ , 

with the same expression for hatted quantities. Then the properties for the g's lead to 

27re 27re 
pi{l/e) ~ _(J, + J2-J3), p2(l/e) ^ -j={Ji + J3 - J2) (6.20) 

27re 27re 
P3(l/e) ~ — (J2 + J3-J1), P4{l/e) — (Ji + J2 + J3) 

and 

#fc+#i=27rna , xGCq. (6.21) 
Furthermore x <-> 1/x relates different sheets of the surface 

^1,2(1/2;) = 27r?n - P2,iix), P3^i{l/x) = -Inm - Pi^six) (6.22) 
Again we state the following relation 



^1,2,3,4(2:) =Pl,2,3,4 



Now we compute pt on Cg. From eq.(6.22) we have equations for ^-cut 

i^i +^^2 = 27rm, 2^3 +^4 = -27rm . (6.23) 

Among the equations (6.21) one has 

"§2 — = 2TTrS^\ X e Ci 
#3 - #2 = 27rn'^^\ X e C2 
#4-#3 = 27rn(3), XGC3 

which give us eq.(6.3). Remaining equations are seen in the Bethe ansatz side through the 
introduction of stacks (see discussion in previous subsection.). Finally from (6.18) and (6.20) 
we again obtain (6.17). 
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6.2 Energy and Momentum 

In this subsection let us observe the perfect matching of both energy and momentum from 
classical and Bethe ansatz calculations. Energy is again defined by 

Repeating the same calculations as for 5*0(4) case we see that once again the result depends on 
the behaviour of pe{z) near points z = ±2, eq.(6.15). In fact everything goes exactly as before 
so that we will arrive at the equivalent of (3.38), namely 

E = 27tMkI + I - E E 4"^ 4"^ - 

y a=l,2,3 p 

Let us recall that this was obtained by calculating the a independent integral for a = — 2 + e. 
Had we calculated it for a = 2 — e and we would have obtained 



E = 2ttMkI --EE ^p^^"^ - 



a=l,2,3 p 

Both results should be equal. This gives us not only the value of the energy 

E = ttM (k^ + kI) , (6.25) 

but also the value of the expression in the parentheses which is nothing but the total momentum 
P obtained as before by summing Bethe equations 

^ = ^ E smh(7r0«) = - E E 4"^ -mL = TrM (4 - kI) . (6.26) 

a a=l,2,3 p 

Since we have already identified M = ^ we observe perfect matching between (6.25,6.26) and 
the classics (2.8). 



7 Conclusions and Prospects 

The 0{2n) cr-models considered in this paper cannot be viewed of course as string sigma models 
in its full quantum version, although the 0(6) cj-model in the classical scaling limit perfectly 
describes the compact bosonic sector of classical superstring on AdS^ x , when the Virasoro 
conditions are imposed. We rather consider these quantum theories as toy models bearing many 
realistic features of a quantum string sigma model and shedding some light on the quantization 
of the full Green-Schwarz-Metsaev-Tseytlin superstring. These models show some important 
features of, yet unknown in full detail, string Bethe equations: the nested nature of these 
equations, following from the symmetry algebra of the target space, the underlying discrete 
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quantum degrees of freedom whose role is played in our case by physical particles and by the 
magnon type excitations on a dynamical inhomogeneous lattice created by particles etc. 

Let us point out that the starting point of our construction, the physical S-matrix, can 
be viewed as describing the anti-ferromagnetic states of a specific inhomogeneous lattice spin 
chain, as pointed out in [42]^^, what might be related to the proposals of [37]. The resulting 
inhomogeneous dynamical spin chain of our paper is however treated in the ferromagnetic way 
in order to restore the classical limit. 

This dynamical lattice becomes regular and rigid in the strong coupling limit of o"-model 
giving the XXX-type spin chain, similar (up to some normalization factors) to the spin chain 
describing the 1-loop N = 1 SYM dilatation operator. This observation might lead to the 
right mechanism reproducing the all-loop SYM Bethe ansatz of the papers [26,27,65], including 
the periodicity with respect to momentum due to the regular lattice. The Hubbard model 
interpretation of the paper [65] should be somehow incorporated in the full quantum string. 

The asymptotic freedom of the quantum 0(2n) models cannot be a part of this full string 
CT- model on AdS^ x S^. It should be rather a finite theory, with zero beta function and without 
any logarithmic divergencies in the world sheet perturbation theory. 

What could this string sigma model look like? We find the construction of quantum string 
states out of the physical particles based on the Yang periodicity equation of eq.(3.1) very 
promising. Of course there can be no massive particles there, and the ^-filling of the vacuum 
described here and in [64] should be based not on the confining box potential but on a subtle 
equilibrium between all roots of the Bethe ansatz equations. 

It is plausible to think that the system of Bethe equations for the full theory is based on 
the psu(2,2|4) algebra and may look like an affine extension of it, as we saw in the bosonic 
cases. Even though the scattering of particles in the critical string theory could be qualitatively 
different from that of massive particles, we still expect that such Bethe ansatz equations may 
serve as a quasi-particle description [73, 74] and provide us with a basis for the Fock-space like 
construction of the quantum Hilbert space. 

Apart from the string theory applications, we think that our method of restoring the classical 
limit of integrable sigma models is interesting in itself and can be applied to many problems of 
two-dimensional field theories, especially for the construction of the quasi-classical approxima- 
tion. It is also useful for proving the validity of more heuristic approaches to such field theories, 
like of the Zamolodchikovs' bootstrap S-matrices. In particular, we consider the coincidence of 
our classical limit in the Bethe ansatz with the finite gap solution of the 0{2n) sigma models 
as a decisive proof of correctness of the bootstrap approach. 

Another interesting possibility based on our observations is the construction of quantization 
schemes of integrable models starting from the algebraic curves of their finite gap solutions 
(quantization of KdV system in [75] provides an interesting example) . To quantize the algebraic 
curve, we have to choose its particular projection (Riemann surface) and introduce the Bethe 
^®We thank D. Serban for the discussion on this issue. 
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roots which form the cuts of the Riemann surface in the classical limit. The quantization might 
be possible in different projections, but some of them should be more convenient than the others. 
We saw for example from our construction that quantization of the finite gap algebraic curve of 
0{2n) sigma model looks more natural after a change of projection by Zhukovsky map. 

An interesting application of our methods might be the calculation of quantum 1/L cor- 
rections in the compact sector of superstring. It was pointed out in [33] that the direct string 
calculations of [11] of this effects can be reproduced from the compact sector only by the use of 
a natural (^-functional regularization, without a direct computation of fermionic contributions. 
Our Bethe ansatz is a good starting point for such calculation. 
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Note Added: 

We thank T. Klose and K. Zarembo for informing us on their forthcoming work [76] where they 
also developed an approach to the quantization of stringy sigma models in various subsectors 
of the superstring theory on AdS^ x S^. Their approach seems to be very different from ours, 
although both give the right classical string limit. Hopefully, both quantizations describe the 
same system and provide an important complementary information on it. 

Appendix A Charged Particles in a Box 

Let us analyse Bethe equations eq.(3.18) with no u's nor v's present. Furthermore we will 
considering a single mode number for ^'s (which in the continuous limit means only one 6 cut). 

Has we explained in section 3, the /icoshMvrx potential is essentially a box when /i = 
is small, see figure 1. Furthermore, provided we are considering 6 roots not to close to the edges 
of the distribution, i.e. the walls of the box, coulomb interaction in the thermodynamical limit, 
where L ~ M — > oo, becomes exact. 

To gain some intuition about this system we will find exact distribution of the system of 
Coulomb particles inside square box. For the sake of generality we consider a box will charged 
walls with charge qM. Then the equilibrium condition is 




X 2 j 
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The solutions of this discrete equation, i.e. the positions of the particles are the zeros of the 
Jacobi polynomials [17]. Indeed, defining Q{x) = Ylji^ ~ 



Q"{xi) 1 X 



AqM- 



2 ' 



3^1- 

which implies that the L degree polynomial R{x) = Q"(x)(4 — x^) — AqM xQ'{x) has the same 
zeros as Q{x) so they must be the same up to a multiplicative constant. Comparing the x^ 
coefficient one finds R{x) = {—L {L — 1) — AqM L) Q{x) so that 

Q"(x)(4 - x^) - AqMxQ'ix) + (L(L - 1) + 4q M L)Q{x) =0. 

Thus 

Q{x) oc Pf-'^"- (I) 

where Pn^{z) are the Jacobi polynomials. Particles will be located at the zeros of this polyno- 
mials. 

Let us analyse the large M, L limit. From the differential equation for Q one has, for the 
resolvent G = 

M dx ~ ^ ^^^+^^4_^2^^^^ \m^^ M M^)a-x^- 
Then in the large L ^ M limit 



2,x-2J(2, + A)^2_4(^ + 4^ 
G{x) = V __ _ (A.l) 

Thus, for small g, the density has 2 peaked maxima near the walls at ±2. It vanishes before 
reaching the wall as a square root. In the limiting case, when g — > 0, the support goes to [-2,2] 
as expected and the density becomes 

1 



p{x) 



which reproduces eq.(3.20) for m = in proper normalization. 

We considered this exactly solvable example to demonstrate that influence of the walls results 
in the inverse square root behavior of the density near the walls. 

Appendix B Relating resolvents 

Equivalence between (3.15) and (B.5) will be estabilished in this section. We define resolvents 
in the x and z planes as 

H^{x)^! J^dy^ G^iz)^! P^^du:, y±H = V.±V^) (B.l) 
Jc±x-y Jc± z-w 2 V / 
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where C± are cuts outside (inside) the unit circle in the x plane. Since 

G±{z)= [ ^M^uj'{y)dy= [ ( + 1 + ] dy (B.2) 

Jc±z-w{y) Jc^ yx{z)-y y ^-yj 

we have 

G± (z) = H± [x± (z)] + H± [x^ {z)] - H± (0) . (B.3) 

Derivation of BAE from Classical Equations 

It is natural, see (B.3), to define 

G,{z) = H+{x+) + H+ix^) - H+{0) , 
-Gu{z) = H^{x+) + H^{x^) - H^{0) , (B.4) 

Ge{z) = _ _ + 2H+{x^) - 2g„(x+) - 2H+{0) . 

Vz^ - 4 

These functions are indeed resolvents, i.e. they have only cuts and behave as 1/z at infinity. 
Indeed from x-^{oo) = oo, a;_(oo) = we see that they vanish for large z. Now, d, (Gu) have 
only V (u) cuts since at the cut [—2,2] one has x+ <-> x_. On the other hand Gq^z) has only 
got 9 cut. 

Then, it is easy to see that (3.15) is satisfied if 

m(x) = + + vrn, H(0) = Inm (B.5) 

X — I X + 1 

Derivation of Classical Equations from BAE 

Now let us assume that eq.(3.15) is satisfied and let us show that H{x) defined by 

Hix,) ^ -0„(^) + ^ - + 2™ (B^7) 

2 yjz^ - 4 

satisfies eq.(B.5). To see this in suffices to notice that i— ^=^=^ — x±—i ~^ x±+i ' -l^oreover since 
x_ = means z = 00 we immediately have (0) = 27rm. Since (fig — Gu — G^ = — 27rm, H[x) is 
continuous on the unit circle and is thus a well defined function of x. Furthermore it is obviously 
zero at infinity. 

Appendix C Bootstrap and Bethe Ansatz 

To make the account of this paper self-contained, we decided to review in this appendix the 
original bootstrap program of [38] for the SO{A) = SU{2) x SU{2) non-linear sigma model. We 
follow closely, apart from some minor changes, the original notation. Furthermore we will use 
the obtained 5'-matrix to carry the algebraic Bethe ansatz procedure and quantize the sigma 
model [71]. 
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Figure 9: Pictorial representation of Yang Baxter equation. In each diagram time flows in the up direction and 
initial particles are aligned, from left to right, according to its rapidity (the faster the more to the left). Consider 
for instance the last diagram in the Ihs of the last equation. It represents the annihilation of two blue particles to 
form any possible pair of green particles. One of the particles of this pair will be reflected back by a red particle 
to meet again its green twin and annihilate to give rise to a pair of black particles. The fastest "in" -particles are 
the blue ones while the slowest "out" -particles are the black ones. 

Particles are associated to the noncommutative symbols Ai{9) where i = 1,...,4 stands 
for the isovector index and the rapidity 6 parametrizes the momentum as p = mo cosh(7r0). 
"In" states are identified with products arranged in the order or decreasing rapidity while "out" 
states should be arranged in the opposite order. The transition from a 2 particle "in"-state to 
an "out" -state (assuming no particle creation, which is indeed a condition for integrability) is 
written as 

A{9)A,i9') = S2{A9)(^~H^0)6^,Yl^Akie')Akie) + Aj{e')A{9) + h~H^0)Ai9')^^^ 

where A9 = 9' — 9. In this expression the first term represents annihilation, the second one is 
the transmission contribution and the third one the reflection term. In general, the transition 
of the "in" -state 

Aiie^) Aji92) Akie^) 

to the "out"-state is not so simple. However, for integrable theories where one has a large 
quantity of conserved charges it happens that the existence of these charges imposes that, in 
every region where the particles are far away from each other, the wave function is a sum of 
plane waves where in each region the set of momenta {pi} is the same. Then for these integrable 
theories we can obtain the "out"~state by applying (C.l) consecutively to the pairs in this 
triplet until the A's are ordered by decreasing rapidities. This is called factorizability. Now, 
consistency requires that the final result should not depend on whether we arrive at region 
9^ < 92 < 9i by scattering first Aj at first with either ^4^ or Ai. By doing it in both ways 
and equating both results one obtains for each choice of isospin indexes some constraint on the 
functions in (C.l). These equations can be obtained in a straightforward fashion and have a nice 
physical interpretation - see caption in figure 9 - all other possible choices of isovector indexes 
(i.e. colors) give rise either to the same constraints or to some trivial identities where the Ihs 
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and the rhs are already identically equal. One obtains: 

h{9) + h{e') = h{9 + e') (C.2) 
h(e)+g{e + 9') = g{9') 

plus another similar, yet bigger, equation corresponding to the third identity if figure 9. The 
solution of the 2 equations yield h{9) and g{9) in terms of two constants A and k. When plugging 
these solutions into the third equation one obtains k(A). At the end of the day one has 

h{9)=g{i\-9) = '-^. (C.3) 

Thus, out of 3 functions we are now left with one function and one constant. From (C.l) one 
reads off the S-matrix 

Sil{9) = S2{9) [Sik^ji g~\9) + 6,,5ik + 5uSjk h^HO)] (C.4) 

Changing i ^ j and the channel s = 4mQ cosh^(7r0/2) u (i.e. 9 ^ i — 9) should leave the 
S-matrix invariant. This is crossing-symmetry. It implies 

S2{B) = S2{i - 9) (C.5) 

and h{9) = g{i — 0), or A = 1. Finally we must impose the most natural requirement, namely, 
unitarity. Setting 

Si^n{-0)ST{0) = 6,j6ki S2{-9)S2{9) (1 + h-\-9)h-\9)) + 6,kSji (...)+ 6uSkj (. . .) 

to be equal to Sij5ki one obtains 3 equations. The exact expressions inside the parentheses are 
not relevant for our discussion. It suffices to say that, for the g and h we found, (C.3), they 
vanish identically. Thus one is left with 

S,{-9)S2{9) = J^. (C.6) 
From (C.5) and (C.6) it follows that S2{9) is given by 



times a CDD factor 

L 

r 

sinh 7x9 — 1 sin 



jj^gl) _ TT sinh 7x9 + i sin ^ 



k=l 



where are arbitrary real numbers. The form (C.7) is needed to have the right pole and zero 
structure according to (C.5) and (C.6) while the ambiguity (C.8) is unfixed since f{9)/f{i — 9) = 
f{9) f{—9) = 1 . Absence of additional bound states forces one not only to exclude this factor 
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but also to introduce the i in Sq{6) [43,77]. For the non-linear sigma model the correct S-matrix 
is indeed given by the minimal choice (C.7). To verify this claim some convincing cross checks 
are done [38,78]. Nevertheless notice that for what concerns the classical analysis we did not 
need to know anything about this CDD factor since for large 0's it is given by 1 + 0{e~'^'^^) 
being therefore irrelevant in the scaling limit. If this CDD factor was present it could, however, 
in principle contribute into quasi-classical corrections. 

Finally since SO{4) = SU{2) x SU{2) we can replace i, k,j, I by (q, d), /?), (a', d'), (/?', /^') 
and write (C.4) as 

{^e e.,e^,e-'^'e-'^' + e^O - ^) 6^ 6^ 6^ - ^{e - ^) €€4^ 



or 



s{e) = SRie) X sue) , SlA^) = So{e)']zie) , n = —— (c.9) 

o — t 

where P is the permutation operator in X C^. Notice that, since 50(4) = SU{2) x SU{2), 
we could have started with this ansatz where 



instead of starting with (C.l). The absence of the channel where two particles annihilate is 
natural since the particles are charged with left/right charge. Consistency relation would be 
represented by the first equation in figure 9 only, or equation (C.2), which states that h{9) is 
linear in 9. Then we would continue by imposing unitarity and crossing-symmetry to arrive at 
the same results we did for the form of 5*0 and the proportionality constant in h{6). The results 
of this appendix could be easily generalized to both 0{N) and SU{N)l x SU{N)r non-linear 
sigma models or to many other groups. 

Let us now start from the obtained S-matrix (C.9) and carry out the algebraic Bethe ansatz 
program [71] . We introduce a ghost particle in the auxiliary space and scatter it through the 
other particles. We want to diagonalize 

L 

Trf{e) , f{9) = l[Sok{9-ek) 

k=l 

where we trace over the auxiliary space. This is a relevant problem because if we solve it for any 
9 and then set 9 = 6k then this means that the ghost particle will change its quantum numbers 
with particle k since Sok{0) = —Pok x Pok- In other words, 

— Trr(6'fc) = Sk,k~i{9k — 9k~i) ■ ■ ■ Sk,i{dk — 9i)Sk,N{9k — 9n) ■ ■ ■ Sk,k+i{9k — 9k+i) 
so that the problem of imposing periodicity of the wave function reads 

_ g*mo£smhMfc) Xrf 1^) = I*) . (C.IO) 
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Let us now perform the diagonalization of C.IO. We consider 

1^) = n BUui) n Br{v3)\ ^{Oi, • • • , Ol)) (C.ll) 
i=l j=l 

where is the state with L particles, where the right and left spin of every particle is pointing 
in the up direction, and 

with a similar definition for the left sector. Acting on VL one has 



1 / 0-§(r^ + l) -§r- \.^. 1 {e-i)\^) 



* 



e + |(r3-l) 6-i\ e\n) 



The upmost right element is not important for our discussion. However the vanishing of the 
element in the left down corner is. It implies that | O) is eigenvalue of both A and D with 
eigenvalues 

L 

A{e)\n) = J{So{e-e^)\n) 

a=l 



OL = \ 



So now we have to understand how A and D pass through the S's. Consistency relations can 
be cast as Sx2(6)Sxz(Q + Q')S2-i{Q') = S23{0')Si3ie + e')Si2{9) which imply (cf. [79]) 

n,Li(^) T^dO') Saa' iO' -e)= Saa'iO' - 6) T^^^O') T^dO) ■ 

where a and a' are two auxiliary spaces. This gives us the commutation relations between 
the elements of the transfer matrix (C.12). In particular 

[B{e),B{e')] = 

A{e)B{e') = ~ ' B{9') A{e) + ^ B{6) A{9') (C.13) 

D{e)Bie') = + ' B{e') D{e) - ^ B{e) D{e') . (c.i4) 

for symbols in the same right or left sector. Symbols in different sectors commute to zero of 
course. Then, acting on (C.ll), one has 

-T,f{e)\^) = {AR{e) + bR{e)) x {AL{e) + bLm^BL{ui) ^\{BR{vj)\vi{ei,...,eL)) 

i=i i=i 

a=l \t=l «=1 a=l / 

\t=l 1=1 a=l / 
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where dots stand for undesirable terms which would make |^) not to be an eigenvector of TrT 
while the displaced term is the one we obtain ignoring the second term in the rhs of both (C.13) 
and (C.14). The condition that these undesirable terms vanish gives us a set of equations for Ui 
and Vj. There is however a shortcut to arrive at these equations provided we know that these 
terms can indeed be killed. The argument is the following - each of the two last terms inside the 
big parentheses came from the diagonalization of a product of Q = 6 — iP. The diagonalization 
of such a product of operators must yield a polynomial in 6 therefore the residues of the apparent 
poles which seem to be part of the eigenvalue for 6 = Ui (or vj) must vanish. This implies 

^ ^ J^ n,-u, + i ^ Uj-Oa ^^.^^^ 
Uj -Ui-l Uj -6a + i 

1 = fr!i^^-rT^i_^ (C.16) 

Furthermore (C.IO) reads 

^.^oZ:sinh.., -Q s2^Q^ _ -Q 0^_U^ ^ e^_V^ ^ ^ ^^^^^^ 

a^P 1=1 ^ 1=1 ^ 

Equations (C.15-C.17) coincide precisely with (3.4-3.6) after the trivial shift {u,v) — > (u — 
i/2, V — i/2) in the former. 

Appendix D Bethe Ansatz Equations for 0(2n) Sigma-Model 

Zamolodchikovs' S-matrix for the 0(2n) sigma-model takes the form [38,80] 



fa'b' _ ^pa'b' ^ T>a'b' 

ab Q ab i(^n-l)-e 



(D.l) 



where the overall factor is given by 



= r(A + ,)r(i-,)r(| + ,)r(^ + A-,) 
^ r(i + A - ^)r(^)r(i - ^)n\ + a + ^) ' ^ ^ 

~ 2n - 2 ' 2n-2' ^^'^^ 

and 

iit=5i5'b. Pab=^aSb, =5ab5''''' (D.4) 

are 0(2n)-invariant tensor bases. Let us restrict ourselves to n > 2. The above form is deter- 
mined as a minimal solution satisfying the Yang-Baxter equation 

the unitarity condition 
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and the crossing symmetry 

Sit{Kn-l)-e)=Sil,{9). (D.7) 

By non-trivial comparison with perturbation theory [38,80-82], this S-matrix is beUeved to 
describe the scattering of vector particles in the 0{2n) sigma-model. Let us consider the system 
of L particles put into a one-dimensional periodic box. The system is described by the Bethe 
wave function ip"''^"'"'^{9i, . . . , 9l) where {Oa} are rapidities associated with each of the particles. 
The wave function reflects the elastic, factorisable property of the scattering as 



It also satisfies the periodic boundary condition 



{9i,92, 9s,. . . , 9l). 



a2---a]^a\ 



{92, . . . ,9l,9i 



{9i,...,9l). 



(D.8) 



(D.9) 



The momentum of the ith particle is given by p{9a) = /isinh(^^0Q,). Following from these 
properties, the wave function obeys the Yang equations 



for a = 1, ... , L. We introduced the transfer matrix T{u)a\l.'^aL, given by the trace 



where 



^"aiut^-tl = Sr,\ (n - 9^)Sllll {u-92)--- Sl^ll (u - 9l) 



^fe b, 



(D.IO) 

(D.ll) 
(D.12) 



is the monodromy matrix. 

The matrix- valued equations (D.IO) can be diagonalized with the help of the nested Bethe 
ansatz [71]. Supplied with an appropriately defined Bethe vector (wave function) [83], the 
transfer matrix T{u) is diagonalized independently of the value of u. The lowest eigenvalue is 
given by 

n 

T{u) = W{u)Y,{tk{u)+tk{u)), (D.13) 



where 



tk{u) 


= Yk- 


-i{u 


tn-l{u) 


= Yn- 


-2{u 


tn{u) 


= Yn- 


-i{u 


tk{u) 


= Yk- 


.i(n 


tn-liu) 


= Yn- 




in{u) 


= Yn- 


-i(m 



u 



k=l 



{k-l)i 



) , for k = 1, . . . ,n 



(n-l)i^-ly^ 
2 J Yn(u - 2 



2 

ni\~l 



(n~2)i\ 



{n~2)i 



i2n-2--k)i ^Yk{u - for k = l,...,n-2. 



(n-l)^ 
2 



^^)Y^{u-W\ 



(D.14) 
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= n(0) ^/(l) 

o^<yo O 



n 



(n-l) 



u 



in) 



Figure 10: Dynkin diagram for the Bi^'^ root system 



Scattering terms W{u),Yk{u) are defined as 



Ko 

Wiu) = lla2{u-u^^^), (D.15) 



Kk (k) I i 



^^(^) = n k k = 0,...,n. (D.16) 

We sometimes use the notations := 9a and Kq := L for the sake of simplicity, {u^''^} are 
rapidities of magnons introduced in the construction of the Bethe vector [83]. 
The diagonalized Yang equations are then given by 

^-.,sinh(^..) ^ -Q _ -Q ° ^^^^ + 2 ^ (D.17) 

P^OL j=l — Uj — 2 

^(B)-^-\(f))- r(A + (^)r(i-c^)r(i + v.)r(i + A-^) 

5o(^) - - -r(A-^)r(i + ^)r(i-v.)r(i + A + ^)- ^^'^'^ 

To make the diagonalization complete, the rapidities {u^'^^} have to satisfy the Bethe ansatz 
equations. Ignoring the common scalar factor W{u), T{u) is essentially a polynomial in u by 
construction and thus has no pole at any finite value of u. The requirement of cancellation of 
poles between two consecutive terms in T(n) gives rise to the Bethe ansatz equations 

XJ-ii (fc) (0 j/ I 

i=0j=iul ' -u)' - ^{ak\ai) 

{afc} {k = 1, . . . ,n) denote the simple roots of the Dn = so(2n) Lie algebra and (a^la/) expresses 
inner product defined in the root space. The non-zero elements read explicitly {ak\o(k) = 2 
(fe = 1, . . . ,n), {ak\ak+i) = — 1 (A; = 1, . . . , n — 2) and (an-2|on) = —1, which are encoded in 
the Dynkin diagram (see Fig. 10). In addition we formally introduced an extra root oq giving 
(ao|ckfc) = — <^ifci so that the interaction between O^s and uf^^^s are included in the common 
notation (D.19). Note that the norm of ao can be fixed as (ao|ao) = 1) if we take account of 
the asymptotics of the self interaction between ^^'s: Observing the asymptotics 



So{e) = i-- + o{e-^), (D.20) 
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we see that the self interaction between 9a s is half the strength of that of Uj's. In the Sutherland 
limit, the Yang equation (D.17) and the Bethe ansatz equations (D.19) are unified into the form 



^iak\ai)(^i{uf'^/M) =27rni^\ for k = 0,...,n, (D.21) 

1=0 



(k) 

with mode numbers n ■ G Z and resolvents 



Note that oq precisely corresponds to the affine root which extends the root system of Dn to 
that of the Bn^ affine Lie algebra. This enhanced symmetry may be viewed as a part of the 
symmetry of the conformal field theory arising at the renormalization fixed point. 
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